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ABSTRACT 

Recent astronomical observations indicate that our Universe is undergoing a 
period of an accelerated expansion. While there are many cosmological models, 
which explain this phenomenon, the main question remains which is the best one 
in the light of available data. We consider ten cosmological models of the accel- 
erating Universe and select the best one using the Bayesian model comparison 
method. We demonstrate that the ACDM model is most favored by the Bayesian 
statistical analysis of the SNIa, CMB, BAO and H(z) data. 



1. Introduction 



Recent observations of type la a supernovae (SNIa) provide the main evidence that the 
current Universe is in an accelerating phase of expansion (IRiess et al. 19981 : iPerlmutter et al. 19991 ). 
There are many different cosmological models used in explanation of an accelerating phase of 
evolution of the current Universe. They can be divided into two groups of models according 
to 'philosophical' assumptions on a cause of the accelerated expansion. In the first type of 
explanation the conception of mysterious dark energy of an unknown form is used, while 
in the second one it is postulated some modification of the Friedmann equation. Here we 
choose five models which belong to the former group as well as five ones which belong to the 
latter one. All the chosen models are assumed to be spatially fiat. 



If we assume the Friedmann-Robertson- Walker (FRW) model in which effects of non- 
homogeneities are neglected, than acceleration can be driven by a dark energy component 
X (matter fiuid violating the strong energy condition). This kind of energy represents 
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roughly 70% of the matter content of the present Universe. The model with the cosmological 
constant (the AC DM model) has the equation of state for dark energy as follows: px = —px 
(IWeinberg 19891 ) . The r nodel with pha r itom dark energy has p y = wxPx, where wx (< — 1) 
is a negative constant (ICaldwell 20021 : iDabrowski et al. 20031 ). The next one is the model 
with a dynamical coef ficient of the equation of state, paramet erized by the scale factor a: 
w{a) = wo + wi{l — a) (IChevallier fc Polarski 200ll : iLinder 20031 ). The other simple approach 
is to represent dark energy in the form of a minimally coupled scalar field (p with the potential 
V{(p). In cosmology the quintessence idea is important in understanding a role of the scalar 
f ield in the current Universe. We consider the power-law parameterized quintessence model 
( jPeebles fc Ratra 19881 : iRatra fc Peebles 19881 ) . In this case density of dark energy changes 



with the scale factor as px = Pxocl where wxio) is the mean of a coefficient of 

the equation of state in the logarithmic scale factor 

/ wx{o)d\n.a 



wx[a 



Jd{lna) 



and has the following form wx = woa"" (iRahvar fc Movahed 20071 ). The first group is com- 
pleted with the model with the generalized Chaplygin gas, where Px = (here v4 > 

Px 

and a = const). We gathered above models together with their Hubble functions (with the 
assumption that the Universe is spatially fiat) in Table [H 



As we have written before we also consider five models offering the explanation of cur- 
rent acceleration of the Universe in an alternative way to dark energy. The brane models 
has postulated that the observer is embedded on the brane in a l arger space in wh ich grav- 
ity can propagate: the D yali-Gabadadz e-Porrati model (DGP) (IDvali et al. 20001 ). Sahni- 
Shtanov brane 1 model (IShtanov 2000l ). The Cardassian model, in which the Universe 
is flat, is matter dominated and accelerating as a consequence of the modiflcation of the 
Friedmann flrst integral as follows 3H^ = P + -Bp", where B is a. constant and the energy 

Table 1. The Hubble function for cosmological models with dark energy 



model 



(2:)relation 



1 ACDM model H'^ = H^^{nm.Qi^ + zf + i'i- - ^1^,0)} 

2 model with generalized Chaplygin gas = + zf + (1 - nm,o)[As + (1 - As){l + 2)3(1+°^)] TT^ | 

3 model with phantom dark energy H'^ = _H'|{r2ni,o(l + + (1 - nm,o)(l + 2)«(i+™^)} 

4 model with dynamical E.Q.S //^ = ^2 jn,^ q(i + + (1 - nm,,))(l + 2)3(i"o+»i+i) exp[-^^]| 

5 quintessence model = Hi |Oni,o(l + ^)^ + (1 ~ ^^m,o)(l + 2)3(i+™o(i+z)"°) | 
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density contains only dust matter and radiation (IFreese fc Lewis 20021 ). We also include 
in the analysis th e bouncing model arising i n the context of loop quantum gravity (the 
BACDM model) Jsingh fc Vandersloot 20051 : ISzvdlowski et al. 2005h and the model with 
energy t ransfer between the da rk matter and dark energy sectors (the A decaying vacuum 
model) (ISzydlowski et al. 20061 ). We gathered above models together with their Hubble 
functions (with the assumption that the Universe is spatially flat) in Table [2l 

The main goal of this paper is to compare all these models in the light of SNIa, CMB, 
BAO and H(z) data. We use the Bayesian model comparison method, which we describe in 



lection (see e.g. Liddle 2004: John & Narhkar 2002: Saini et al. 2004 


: Parkinson et al. 2005: 


Mukheriee et al. 2005: Beltran et al. 20051: Mukheriee et al. 2005a: ^ 


>zvdlowski & God 


owski 2006: 


Godlowski & Szvc 


lowski 2OO5I: Szvd: 


owski et al. 2006a: lliddle et al. 20061: Liddle 2007 


Sahlen et al. 2007 


Serra et al. 200?!: 


Kunz et al. 20061: 


Trotta 2OO7I: Trotta 2007a). Recentlv the Bavesian in- 



formation criteri a were applied in t he context of choosing an adequate model of acceleration 
of the Universe (jPavis et al. 20071 ). The authors showed preference for models beyond the 
standard FRW cosmology (so-called exotic cosmological models) whose best fit parameters 
reduce them to the cosmological constant model. 



2. Model comparison in Bayes theory 

Let us consider the set of K models: {Mi, ■ ■ ■ , Mk}- In the Bayes theory the best model 
from the set under consideration is this one which has the largest value of the probability in 
the light of the data [D), so called posterior probability 

, , , P(D\Mi)P(MA 
PmO) = ^ ^p^'^^^ (1) 

P{Mi) is the prior probability for the model indexed by i, which value depends on our 
previous knowledge about model under consideration, that is to say without information 
coming from data D, and P{D) is the normalization constant. If we have no foundation to 
favor one model over another one from the set we usually assume the same values of this 
quantity for all of them, i.e. P(Mj) = -^, i = 1, ■ ■ ■ , K. 

To obtain the form of P{D) it is required that a sum of the posterior probabilities for 
all models from the set is equal one 

K K 
1=1 i=l 
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Therefore conclusions based on the values of posterior probabilities strongly depend on the 
set of models and can change when the set of models is different. 

P[D\Mi) is the marginal likelihood (also called the evidence) and has the following form 



P(D\Mi 



L{ei\D,Mi)p{ei\Mi)dei = Ei, 



(2) 



where L{9i\D, Mi) is the likelihood of the model under consideration, 9i is the vector of the 
model parameters and P{9i\Mi) is the prior probability for the model parameters. 

In the case under consideration we cannot obtain the value of the evidence by analytical 
computation. We need a numerical method or an approximation to this quantity. 



Schwarz (author?) ( ISchwarz 19781 ) showed that for iid observations {D = {xi},i 
1, - ■ ■ , N) coming from a linear exponential family distribution, defined as 



f{xi\9) = exp 



Y,MdMx,)+b{9) 



,k=l 



S = d, 



where wi, . . . ,ws,b are functions of only 9 G R'^, ti, . . . ,ts are function of only Xi, the 
asymptotic approximation (A^ —* oo) to the logarithm of the evidence is given by 



ln£; = ln/:-^lnA^ + 0(l), 



(3) 



where C is the maximum likelihood and 0(1) is the term of order unity in A^. In this case 
the likelihood function has the following form 



L{9\D,M)=Ulj{xi\9)=exp 



Nlj2Md)Tk{D) + bi9) 



,k=l 



where Tk{D) = jj: J2iLi'^k{xi)- The integral ([2]) can be writing as 



exp[Ng{9)]P{9\M)d9, 



(4) 



where g{9) = J2k=i''^k{0)Tk{D) + b{9). This integral has the form of the so called Laplace 
integral. Assume that g{9) has maximum in and P{9q\M) ^ 0. When N oo exp[Ng{9)] 
will be a sharp function picked at ^o- Then the main contribution to integral (jlj) comes from 
the small neighborhood of 9q. In this region P{9\M) ^ P{9o\M), we can also replace g{9) 
function its Taylor expansion around 9o: g{9) = g{9o) — ^{9—9o)^C^^(9—9o), where fC~^l,„. = 



and extend the integration region to whole R . One can gets the asymptotic 
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of the integral E = exp[Ng{eo)]P{9o\M) {^Y VdetC and In^ = Ng{9o) -jlnN + R, 
where R is the term which not depend on N. One can see that Ng{6o) = In L{6o\D, M), 
where Oq is the point which maximize g{9) = \n L{9\D , M) , so is equivalent to 9mle (the 
maximum likelihood estimator of 9). Finally one can obtain result ([3]). 

According to this result Schwarz introduced a criterion for the model selection: the best 
model is that which minimizes the BIC quantity, defined as 

BIC = -2\nC + d\nN. (5) 



This criterion can be derived in such a way that it is not required to assume any specific 
form for the likelihood function but it is only necessary that the likelihood function satisfies 
some non-restrictive regularity conditions. Moreover data do not need to be independent and 
identically distributed. This derivation requires to assume that a prior for model parameters 
is not equal to zero in the neighborhood of the point where the likelihood function under a 
given model re aches a maximum and tha t it is bound in the whole parameter space under 



consideration (jCavanaugh fc Neath 19991 ). It should be pointed out that an asymptotic 
assumption is satisfied when a sample size used in analysis is large with respect to the 
number of unknown model parameters. 

It is useful to choose one model from our models set (here indexed by s) and compare 
the rest models with this one. We can define ABICis quantity, which is the difference of the 
BIC quantity for the models indexed by i and s: ABICis = BICi — BlCg and present the 
posterior probability in the following form 



exp{-^ABIC,,)P{Mi 
Ef=iexp(-lAS/C,,)P(Mfc) 



pmD) = z... - (6) 



Let us assume that we have computed the probabilities in the light of data D for models 
from the set under consideration. Then we gathered new data Di and want to update the 
probabilities which we already have. We can compute probabilities in the light of new data 
using information coming from previous analysis, which allow us to favor one model over 
another: we can use posterior probabilities for models obtained in earlier computations as a 
prior probabilities for models in next analysis. 



We apply this method in evaluation the posterior probabilities for models described in 
the previous section using the information coming from SNIa, CMB, BAO and II{z) data. 
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3. Application to cosmological models comparison 



We start wit h the N = 192 sample of SNIa flRiess et al. 20071 : IWood-Vasev et al. 2007 
Davis et al. 20071 ). In this case the likelihood function has the following form 



L oc exp 



It 



theor 
i 



obs\2 



i=l 



where CTj is known, ^"'^^ = rrii — M (mj-the apparent magnitude, M-the absolute magnitude 
of SNIa), nf"""^ = 5 logio ^Li + M, M = -blogioHo + 25 and Dn = Hodu, where du is the 
luminosity distance, which with assumption that /c = is given by 



d 



Li 



'1 + Zi)c 



dz' 
Hiz' 



In this case we used the BIG quantity as an approximation to the minus twice logarithm 
of evidence and assumed that all models have equal values of prior probabilities. Based on 
our previous experiences we have used following assumptions for models parameters values: 
Hq E (60, 80) for all models and additional: 

• model 2: As G (0,1), a G (0,1) 

• model 3: Wx G (—4, —1) 

• model 4: wq G (—3,3), Wi G (—3,3) 

• model 5: Wq G (—3,3), a G (0,2) 

• model 7: Qn,o G (0, 1), n G (3, 10) 

• model 8: Qint,o G (-1, 1), n G (-10, 10) 

• model 9: n G (-10, 10) 

• model 10: Vlifi G (0,4), ilAbfl G (-1,4) 



We separately consider cases with Qm,o G (0,1) and fim.o G (0.25,0.31). We treat the Hq 
parameter as a nuisance parameter, i.e. we marginalized the likelihood function over this 
parameter in the range assumed before. Posterior probabilities are obtained using equation 
[6l We analyse three sets of models: 1. set of models with dark energy (Tabled]), 2. set of 
models with modified theory of gravity (Table [2]), 3. set of all models (Tabled] and Table [2] 
together) . 
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The results for the case with Qra,o G (0,1) are presented in Table 3, Table 4 and Table 
5 for set 1, set 2 and set 3 respectively. One can conclude that in the light of SNIa data the 
ACDM is the best model from the set of models with dark energy as well as the best one 
from the all models under consideration, the DGP model is the best one from the group of 
models with modified gravity. 

The results for case with Qrn.o ^ (0.25,0.31) are gathered in Table 6, Table 7, Table 8 
for set 1, set 2 and set 3 respectively. The conclusion changed for the set of models with 
modified gravity: here the best one is the Cardassian model. 

In the next step we included information coming from CMB data. Here the likelihood 
function has the following form 



L oc exp 



^^theor ^obs'\2 



where R is s o called shift parameter, R^^''°' = g"^^ w r)^^^ and R°^'' = 1.70 ±0.03 for 



z^ec = 1089 (jSpergel et al. 20061 : IWang fc Mukherjee 20061 ). It should be pointed out that 



the parameter R is independent of Hq. 

The values of the evidence were obtained by the numerical integration. We assumed fiat 
prior for all model parameters. It is known that evidence depends on the prior probabilities 
for model parameters. Assumptions for the model parameters intervals, which we made in 
previous analysis could be not appropriate here. Due to this we made a stricter analysis 
for models with parameters which interval width exceeds one. This width was used for 
convenience. We computed the evidence for different parameter intervals, which do not 
exceed the range assumed before and with a minimal width equal to one. There are of 
course extremely many possibilities. We limited our analysis to intervals (a, b), where a and 
b are integer. Finally we chose the case with the greatest evidence. 

We consider the situation with Qm,o G (0,1) and Qrn,o G (0.25,0.31). The range for 
parameters which change after stricter analysis in the first case: 

• model 3: Wx € (—2, —1) 

• model 4: Wq G (-1,0), Wi G (-2,0) 

• model 5: Wq G (-3,-2), a G (1,2) 

• model 7: ^^,0 G (0, 1), n G (3,4) 

• model 8: fiint,o G (-1,0), n G (-10,-9) 
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• model 9: n G (0, 1) 

• model 10: ^]^,o e (0, 1), QAb^o G (0, 1) 

and in the second case: 



• model 3: Wx G (—2, —1) 

• model 4: wq G (-1, 0), Wi e (-2, -1) 

• model 5: Wq G (-2, -1), a G (0, 1) 

• model 7: fi„,o e (0, 1), n G (3,4) 

• model 8: fiint,o e (-1,0), n G (-2,-1) 

• model 9: n G (0, 1) 

• model 10: Qi^o G (0,1), QAbfi e (3,4) 

Posterior probabilities were obtained using equation [H Here we treated posterior prob- 
abilities evaluated in analysis with SNIa data as a prior probabilities. Results are gathered 
in tables like in previous analysis. We also show the values of the posterior probabilities 
obtained for the intervals assumed at the beginning (numbers in the brackets). As we can 
see the ACDM model is still the best one from the models with dark energy (for both ranges 
of l^m.o)- The conclusion is the same for the set of models with modified gravity: the DGP 
model is the best one in the first case and the Cardassian model in the second. When we as- 
sume that fim.o £ (0, 1) there is no evidence to favor the ACDM model over the DGP model 
(they have the same values of the posterior probabilities) while when we restrict fim.o range 
to (0.25, 0.31) the ACDM model still stays as the best one, with even greater probability. 

As the third observational data we used the meas urement of the baryon acoustic oscil- 
lations (BAO) from the SDSS luminous red galaxies (lEisenstein et al. 20051 ). In this case 
the likelihood function has the following form 



L (X exp 



^^theor ^obs^S 



2al 



where A 



theor 



m,0 



H{z) 
Ho 



J_ f^A Hp 

ZA Jo 



and A""' = 0.469 ± 0.017 for za = 0.35. 



^ZA Jo H(z)_ 

Here values of the evidence are obtained by the numerical integration. We made the 
analogous analysis with the parameters intervals as above with the additional requirement: 
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obtained intervals must at least cover the intervals obtained in previous analysis. For most 
of models the conclusions are the same as in the previous analysis. Below we wrote the cases 
where the intervals have changed 
for the case with flrn,o ^ (0, 1): 

• model 5: Wq G (—3, 1), a G (1, 2) 
and with Q^^q G (0.25,0.31): 

• model 8: Qi^tfi G (-1,0), n G (-2, 1) 

Here we used posterior probabilities obtained in analysis with the CMB data as prior 
probabilities. Results were again presented in described above tables. The conclusion is 
different for the set of all models in the case with ^2^,0 £ (0, 1): the DGP model becomes 
the best one from them. 

Finally we use d the obse rvational H(z) data (A^ = 9) from Simon et al. (author?) 
(ISimon et al. 20051 ) (see also ISamushia fc Ratra 20061 : IWei fc Zhang 20071 . and references 



therein). These data based on the differential ages (^) of the passively evolving galaxies 
which allow to estimate the relation H(z) = - = — t^^- Here the likelihood function has 

\ ' a 1+2 at 

the following form 



L oc exp 



i=l 



at 



where H{z) is the Hubble function. Hi, zi are observational data. 

The values of the evidence were obtained by the numerical integration. As above we 
assumed fiat prior probabilities for models parameters. The ranges for them which changed 
after analogous to previous analysis are as follows: 
for case with fim.o ^ (0, 1) and G (60,80): 

• model 5: G (-3, 1), a G (0, 2) 

• model 8: fijnt.o e (-1,0), n G (-10,0) 

and for case with Vl^^q G (0.25,0.31) and G (60,80): 



model 4: G (-1,0), Wi G (-2,0) 
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Values of posterior probabilities obtained in analysis with the BAO data were treated 
as prior probabilities it this analysis. The results are presented in tables. 

As one can see in the case with Q^.o £ (0.25, 0.31) the ACDM model is the best one 
from the set of models with dark energy as well as the best one from all models considered 
in this paper. The conclusion is different in the set of the models with modified gravity: 

after the analysis with observational H{z) data the DGP model becomes the best one. In 
the case with ri,n,o G (0, 1) the ACDM is still the best model from the set of models with 
dark energy while the DGP model is the best one from the set of models with the modified 
theory of gravity as well as the best one from all models considered. 

As one can conclude final results coming from computation including the shrinking 
parameters interval procedure are the same as final results coming from computation with 
parameter intervals assumed at the beginning: the best model (in each set) does not change, 
but the probability of being the best one are greater for second case. 

As we have written before we used the BIG quantity as an approximation to the minus 
twice logarithm of the evidence for the SNIa data. This approximation gives good results if 
in the set under consideration is one favoured model. The problem appears when we have 
two favoured models with nearly the same probabilities. Such a situation is in the set of all 
models with Qm,o 0, 1 >. ACDM model and DGP models have nearly the same values of 
model probabilities which are the greatest ones in considered set. This enforce us to compute 
the evidence by the numerical integration for these two models (for SNIa data) to check if 
our previous conclusion is true. In Table 9 we gathered the values of probabilities obtained 
after computation of the full Bayesian evidence (case 1) as well as such values obtained when 
the BIG approximation was used (case 2). 

As one can see the conclusion changed. The AGDM model is better than the DGP 
model in the light of all data sets used in paper. The BIG approximation is not enough in 
this cases, gives us wrong answer. 

Finally we compare considered sets of models treating all described before data sets as 
N — 192 + 1 + 1 + 9 independent data. In this case the hkelihood function has the following 
form 



L oc exp 



(/?"""" - 7?"^")' (A"""" - .4"'")' ^ (//(-,) - H,iz,yf 




Here we assumed fiat prior for model parameters in the range described at the beginning of 
this section. We used the BIC as an approximation to —2lnE. The results for the cases 
with Qm.o 0, 1 > and Qni,o 0.25,0.31 > for all described before sets of models are 
gathered in Tables 3, 4, 5, 6, 7, 8. 
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Table 2. The Hubble function for cosmological models with modified theory of gravity 



model 



(2)relation 



DGP model 

BACDM model 
interacting model with A 



^2 = Hg|[v/n;;^;;^(rTiFTn;^+ v'n;^]^! 

_ (i-n„,o)^ 

H'^ = Hi {f7m,o(l + z f - n„,„(l + 2)" + 1 - Cm,,) + f^„,o} 

H2 = H2{Q„,,o(1 + zf + ni„t,0(l + 2)" + 1 - f^m.O - f^int.o} 



H^ = Hl{ f2,,o(l + + nm,o(l + zf 



TTI + (1 + ^) 



-4+4n / i-nr,o-n„,o \ I 



9 Cardassian model 
Q^O = 10-" 

10 Sahni-Shtanov brane I model H'^ = |Qm,o(l + 2)^ + f2<T,o + 2f2;_o - 2^Qi,o V^m,o(l + ^)'^ + f^<T,o + ^1,0 + f^Ab.o} 



,0 



Table 3. Posterior probabilities for models from Tabled] Q^ifl G (0, 1) 



model prior posterior SNIa posterior + CMB posterior + BAO posterior + H(z) posterior SNIa+CMB+BAO+H(z) 



1 


0.20 


0.91 


0.92 (0.95) 


0.91 (0.97) 


0.94 (0.99) 


0.84 


2 


0.20 


0.01 


0.01 (0.01) 


0.01 (0.01) 


0.00 (0.00) 


0.02 


3 


0.20 


0.07 


0.06 (0.04) 


0.07 (0.03) 


0.05 (0.01) 


0.06 


4 


0.20 


0.01 


0.01 (0.00) 


0.01 (0.00) 


0.01 (0.00) 


0.04 


5 


0.20 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.04 



Table 4. Posterior probabilities for models from Table [21 Qui.o ^ (0, 1) 



model prior posterior SNIa posterior + CMB posterior + BAO posterior + H(z) posterior SNIa+CMB+BAO+H(z) 



6 


0.20 


0.89 


0.92 (0.97) 


0.92 (0.98) 


0.93 (0.98) 


0.07 


7 


0.20 


0.01 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.03 


8 


0.20 


0.01 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 


0.13 


9 


0.20 


0.08 


0.07 (0.01) 


0.07 (0.00) 


0.06 (0.00) 


0.74 


10 


0.20 


0.01 


0.00 (0.01) 


0.00 (0.01) 


0.00 (0.01) 


0.03 
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Table 5. Posterior probabilities for models from Tabled] and Table [21 f2m,o ^ (0, 1) 



model 


prior 


posterior SNIa 


posterior + CMB 


posterior + BAO 


posterior + H(z) 


posterior SNIa+CMB+BAO+H(z) 


1 


0.10 


0.51 


0.46 (0.48) 


0.44 (0.46) 


0.43 (0.45) 


0.74 


2 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.02 


3 


0.10 


0.04 


0.03 (0.02) 


0.03 (0.01) 


0.02 (0.00) 


0.05 


4 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.04 


5 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.03 


6 


0.10 


0.39 


0.46 (0.48) 


0.47 (0.52) 


0.50 (0.54) 


0.01 


7 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.005 


8 


0.10 


0.01 


0.01 (0.01) 


0.02 (0.01) 


0.02 (0.01) 


0.01 


9 


0.10 


0.04 


0.04 (0.01) 


0.04 (0.00) 


0.03 (0.00) 


0.09 


10 


0.10 


0.01 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.005 



Table 6. Posterior probabilities for models from Table [T} Qm.o ^ (0.25,0.31) 



model prior posterior SNIa posterior + CMB posterior + BAO posterior + H(z) posterior SNIa+CMB+BAO+H(z) 



1 


0.20 


0.91 


0.98 (0.88) 


0.99 (0.99) 


0.99 (1.00) 


0.84 


2 


0.20 


0.01 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.02 


3 


0.20 


0.06 


0.01 (0.11) 


0.00 (0.01) 


0.00 (0.00) 


0.06 


4 


0.20 


0.01 


0.01 (0.00) 


0.01 (0.00) 


0.01 (0.00) 


0.04 


5 


0.20 


0.01 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.04 



Table 7. Posterior probabilities for models from Table [2l fim.o ^ (0.25,0.31) 



model 


prior 


posterior SNIa 


posterior + CMB 


posterior + BAO 


posterior + H(z) 


posterior SNIa+CMB+BAO+H(z) 


6 


0.20 


0.19 


0.27 (0.26) 


0.35 (0.50) 


0.45 (0.91) 


0.07 


7 


0.20 


0.04 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.03 


8 


0.20 


0.05 


0.13 (0.04) 


0.24 (0.06) 


0.26 (0.04) 


0.13 


9 


0.20 


0.68 


0.60 (0.70) 


0.41 (0.44) 


0.29 (0.05) 


0.74 


10 


0.20 


0.04 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.03 
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The results in this analysis confirm that the ACDM model is the best one in the set of 
models with dark energy as well as the best one in the set of all models for both ranges in 
flmfl- The conclusion that the Cardassian model is the best one in the set of models with 
modified gravity (for both ranges in fim,o) is with contrary with previous inference where 
the DGP was the best one. The reason of this disagreement is related with the different 
constrains on the flm,o parameter for various data sets for considered model. In Table 10 we 
presented the results of parameter estimation performed for all data sets independently as 
well as for all data sets applied simultaneously for the ACDM, DGP and Cardassian models. 



4. Conclusions 

In this paper we gathered ten models of the accelerating Universe. The five of them 
explain the accelerated phase of the Universe in the term of dark energy while the other 
five explain this phenomenon by the modification of the theory of gravity. We used the 
Bayesian model comparison method to select the best one in the set of models with dark 
energy, in the set of models with modified theory of gravity as well as the best one of all 
of them. The selection based on the SNIa, CMB, BAO and observational H(z) data-we 
treat posterior probabilities obtained in one analysis as prior probabilities in the next one: 
information coming from the previous analysis allow us to favor one model over another. 
We used approximation proposed by Schwarz to the minus twice logarithm of evidence in 
the case with SNIa data, and numerical integration of the likelihood function within an 
allowed parameter space (we assumed fiat prior probabilities for model parameters) in the 
other cases. We consider separately cases with Qrn,o ^ (0,1) and with Qrn,o ^ (0.25,0.31). 
We made a stricter analysis for models with parameters which intervals width exceed one: 
we evaluated the evidence for these models for different parameters intervals with minimally 



Table 8. Posterior probabilities for models from Tableland Table [2l fim,o ^ (0.25,0.31) 



model 


prior 


posterior SNIa 


posterior + CMB 


posterior + BAO 


posterior + H(z) 


posterior SNIa+CMB+BAO+H(z) 


1 


0.10 


0.81 


0.91 (0.82) 


0.96 (0.96) 


0.96 (0.97) 


0.74 


2 


0.10 


0.01 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.02 


3 


0.10 


0.07 


0.01 (0.12) 


0.00 (0.01) 


0.00 (0.00) 


0.05 


4 


0.10 


0.01 


0.01 (0.00) 


0.01 (0.00) 


0.01 (0.00) 


0.04 


5 


0.10 


0.01 


0.01 (0.00) 


0.01 (0.00) 


0.01 (0.00) 


0.03 


6 


0.10 


0.02 


0.02 (0.02) 


0.01 (0.02) 


0.01 (0.02) 


0.01 


7 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.005 


8 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.01 


9 


0.10 


0.07 


0.04 (0.04) 


0.01 (0.01) 


0.01 (0.01) 


0.09 


10 


0.10 


0.00 


0.00 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


0.05 
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width equal one, which do not exceed intervals assumed in the analysis with SNIa data and 
finally chose the best one from them (with the greatest evidence) to the next analysis. We 
compare such results with the results obtained in calculation where we treat all data sets as 
= 192 + 1 + 1 + 9 independent data and use BIG as an approximation to the minus twice 
logarithm of the evidence. 

We can conclude that for the case with G (0.25, 0.31) as well as for case with 

^^m,o e (0, 1) 



• the ACDM model is the best one from the set of models with dark energy as well as 
the best one from the set of all models considered in this paper; 

• the Cardassian model is the best one from the set with models with modified theory 
of gravity. 



The paper was supported by the Marie Curie Actions Transfer of Knowledge project 
COCOS (contract MTKD-CT-2004-517186). The authors are very grateful to the referee 
and to Dr W. Godlowski for helpful discussion and suggestions. 
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